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Abstract 

We show that, in addition to 50(4), the Hubbard model at half filling on a 
bipartite lattice has a group of discrete symmetries and transformations. A 



'^ ' unique Hubbard-Stratonovich decomposition of the interaction term, incor- 

Q ' porating both spin and pseudospin variables on an equal footing, is found in 

which these symmetries are manifestly present. The consequences of this at 

>' 
00 ' the mean field and one loop renormalisation group levels are discussed. 
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The Hubbard model is a simple, generic model of correlated electrons on a lattice. It is 
defined by the Hamiltonian: 



H=Ho + Hi-fi'Y.Mr), 

r,(j 

Hq= - Y^ t„' {cl{r)ca{r') + h.c. 

(rr'),cr 

Hi=UY, n^ij)-- n^ir)-- 



(1) 
(2) 

(3) 



where t^r' connects only nearest-neighbour lattice sites and yu' = /x — y controls the number 
density of the system. U is the interaction energy. 

The model was introduced by Hubbard |l[] as a model of itinerant ferromagnets. Since 
then the repulsive case {U > 0) has been used to describe both ferro- and antiferro- 
magnetism 0,|[ and d-wave superconductivity @J^. For an attractive interaction (U < 0) 
it has been applied to the study of BCS s-wave superconductivity and the transition to 
the Bose condensation of strongly bound pairs [j^. Interest in the Hubbard model on a 2D 
square lattice has been particularly intense since the suggestion by Anderson [Q that it is 
relevant to the description of the high-T^ compounds. 

Despite this interest the Hubbard model remains unsolved except in ID. However, in 
recent years, there have been a number of exact results proved. Yang and Zhang [^ have 
shown, on a bipartite lattice at half-filling and zero field, the existence of an 5*0(4) = 
SUs{2) (g) SUr{2)/Z2 symmetry. The group's action generates rotations of the spin and 
'pseudospin' vectors, S and R: 



S = -^Ur)a^.sM, ^sir) 



R=oHW^^i?W' ^/?W 



C|(r) ' 



e*^-''c}(r) 



(4) 



(5) 



where a is the vector of Pauli matrices and vr is the vector (vr,??, ...) in d dimensions. The 
pseudospin vector comprises superconducting pair creation/annihilation operators and the 
number operator. A consequence of pseudospin rotational invariance is the existence, at 



half filling, of a degenerate superconducting/charge density wave (S-CDW) ordering |TO| , p . 
Ostlund |T2| has given a thorough account of the 5*0(4) symmetry, including an enumeration 
of the possible mean field phases of the model. 

Recently Zhang [|1^] proposed an approximate 5*0(5) symmetry connecting antiferromag- 
netism to (i-wave superconductivity. The proposal offers the possibility of understanding the 
phase diagram of the cuprates and has attracted much attention. The symmetry is not exact, 
unlike 5*0(4), but becomes valid under the action of the renormalisation group. 

In this work we examine exact symmetries. We show that in addition to the continuous 
50(4) there are non-trivial discrete symmetries. These play a similar role to CPT invariance 
in relativistic quantum mechanics. In fact, this analogy is very close, as will be clear below, 
because 5*0(4) is the Euclidean analogue of the Lorentz group, 5*0(3, 1). We show that 
these symmetries, augmented by the addition of the Lieb-Mattis transformation, are very 
useful in resolving an ambiguity inherent in any field theoretic treatment based upon the 



Hubbard-Stratonovich transformation [p^,p^ 



eh^'^Jv[b]e-'^'"-^\ (6) 

where JF is a bilinear Fermi operator and 6 is a Bose field. The use of this identity converts 
the problem of interacting fermions into one of free fermions interacting with bosons. To 
employ this technique the quartic term. Hi, must be written as the square of a bilinear 
operator. The method is ambiguous because there are many ways (actually a continuous 



infinity) of doing this |T6| : 



hi=j: 



^(l-«)S-S + ^(l + a)R-R 



(7) 



Whilst this is formally exact for all a, each choice behaves differently under approximation 
schemes (e.g. mean field theory, RPA, renormalisation group). It is customary to pick a 
decomposition which is most convenient for the problem at hand; for instance a = — 1 for 
magnetism or a = +1 for superconductivity. We will show, however, using the discrete 
transformation properties of the Hamiltonian, that there is a single 'best' decomposition. 



Uniquely, approximation techniques applied to this decomposition preserve all the symmetry 
and transformation properties of the full Hamiltonian. 
It is convenient to introduce a four spinor: 



^(r) = (c|(r),q(r),e*"-"c{(r),e*"-V{(r) 



(8) 



in terms of which both S and R may be constructed. Equivalently a 2 x 2 matrix form is 
possible fl^ . 

We consider three discrete symmetries; charge conjugation, C, parity, V, and time rever- 
sal, T, also the Lieb-Mattis canonical transformation, Z. 

Charge Conjugation, C. At half filling {fi' = 0) on a bipartite lattice, for every state of 
energy e(k) below the Fermi energy there is a corresponding state of energy e(k+7r) = — e(k) 
above it. Under charge conjugation a hole of momentum k is equivalent to a particle of 
momentum k + tt. Interpreting filled states as the Dirac sea and holes as positrons, the 
situation is analogous to the relativistic case. With respect to the four spinor, \E'(r), charge 
conjugation is defined by: 



C 



0-1^ 



10 
0-100 
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Y, 



(9) 
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where Y is the antilinear complex conjugation operator. Its effects on the spin and pseu- 
dospin operators are as follows: 



CR,C — — R, 



(10) 

(11) 

(12) 



Parity, V, leaves S, R and the Hamiltonian unchanged and is of little use to us here. 
Time Reversal, T , reverses the signs of the linear and angular momenta: 



fpf-^ = -p, fsf- 



-s. 



(13) 



Cooper pairs are time reversal invariant as is charge, thus: 



TR 1 ^ it , T R/T = Rz 



(14) 



Using these conditions, we may deduce the form of the time reversal operator acting on 
\E'(r): 
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Y. 



(15) 



The antilinear operator, F, is necessary to ensure that the commutation relations, PaQc 
QaPa = ift- are preserved under T. 

The Lieb-Mattis canonical transformation, Z, is defined by: c^{r) \-^ c^{i^), Cj(r) i- 
e*'^'"c|(r). It acts on the four-spinor as a matrix: 

10 0^ 



xpfr) 



1 
10 
10 



^(r). 



(16) 



It maps spin into pseudospin and vice versa, leaves the non-interacting Hamiltonian, Hq, 
unchanged but reverses the sign of the interaction term, Hj: 



S4^R, 



Hq + Hi <^=^ Ho — Hi. 



(17) 
(18) 



The action of Z on the operators C and T is noteworthy: 



C = ZTZ, 



(19) 



implying that the behaviour of S under T is the same as R under C. 

The Discrete Group of Transformations. Using the matrix forms of C, T and Z we see 
that the set, Q = {±1, ±C, ±T, ±Z, ±CT, ±ZC, ±ZT, ±CTZ}, forms a discrete group 
under multiphcation. It comprises 16 elements, falling into ten conjugacy classes. The 
character table is shown in Table |I[ 

Each of the operators Hq, Hj, S and R may be classified as transforming according to a 
certain representation: Hq transforms as the identity representation, T^^\ Hj transforms as 
r(2) whilst S and R transform as couplets; (J), (J) as T^^) and (^'') as T^^) © T^"^). 

The space of bilinear operators is spanned by S and R. Therefore the quartic operator. 
Hi, can be decomposed into products of S and R. The representations of all possible 
products are generated by: 

(r(3) © r(^) © 2r(5)) © (r(=^) © r(^) © 2r(^)) (20) 

= 6rw©6r(2)©4r(3)©4r(^)©8r(^). 

Discarding all but those products which transform as F^^^ and imposing explicit spin and 
pseudospin invariance we arrive at the unique decomposition: 



hi=j: 



-S-S + -R-R 
3 3 



(21) 



We have constructed this decomposition (corresponding to a = in equation ^) in such 
a way that all its symmetries are transparent and are inherited from the transformation 
properties of S and R. Because of this it enjoys a great advantage over the other choices of 
a. This becomes clear when we consider approximation schemes. 

Mean Field Theory. Using equation |^ we construct a mean field theory for each a: 

U 



^MF ^ ^ 



g(l-a){2S-(S)-(S)-(S)} 



+ |(l + a){2R-(R)-(R)-(R)} 



(22) 



The expectations, (...), are taken self consistently with respect to the mean field Hamiltonian. 
For all a and for U > 0, the solutions are either paramagnetic (T > T^) or antiferromagnetic 

6 



(T < Tc) but in both cases (R) = 0. Similarly for U < 0, there is a S-CDW phase for T < Tc 
and (S) = 0. 

The full theory is invariant under the combined action of Z and U ^—>- —U, therefore we 
expect a good mean field theory to share this property. For instance, the values of (S) (for 
U > 0) and (R) (for U < 0) should not change. We find that this is true for the a = 
theory but not for any other. This can be simply understood by considering the action of Z 
and U 1-^ —U on equation ^. The effective couplings which the fields 'see', j {1 — a) for S 
and ^ (1 + a) for R, are exchanged, the effect of which is to alter the expectation values (S) 
and (R). This has dramatic consequences when a = ±1 in which cases the fields of interest 
acquire zero couplings and drop out of the Hamiltonian altogether. The a = theory avoids 
these problems by having the same coupling, j, for both S and R. 

It is desirable that Hartree-Fock theory be reproduced by the scalar version of equation 

Hi = T. [-U (1 - ") ^' + f/ (1 + a) Rl\ , (23) 

r 

where Sz = |(^| — n^) and Rz = |(ri| + n^ — 1). Hartree-Fock is only recovered for the 
a = decomposition. For a ^ there are extra terms, n^j {ua)-, which violate the Pauli 
exclusion principle. These arise because mean field theory does not respect idempotency, i.e. 
rtfjUfj = rifj but n„ {n„) ^ n„. The a = theory was derived without the use of idempotency 
and so reproduces Hartree-Fock theory automatically. 

Landau- Ginzburg Theory. Because the a = decomposition retains both S and R it 
has a further advantage in that, combined with the use of \l/(r), it allows us to study the 
renormalising effects of magnetism on superconductivity and vice versa. To see these we 
construct a Landau- Ginzburg (LG) functional. 

Using equation (|21| ) in the Hubbard Stratonovich transformation and then integrating 
out the fermionic degrees of freedom leaves an effective bosonic action in terms of the fields 
$ and A which couple to S and R. In the static limit the quadratic part reads: 



s^'^ = lT.(^-Txi<i))^i<i)-H-q) 



(24) 



^E(l + fx(g))A(g)-A(-g), 



where the susceptibihty is: 



xiQ) 



E 

k 



/(Ck+q) - /(Ck) 



(25) 



Ck+q — Ck 

To construct the LG functional we need only long wavelength (small q) excitations about 
stable modes. The sign difference in the expressions for $ and A means that (for U > 
0) the stable modes are at vr (antiferromagnetism) and (r^-pairing [l^) for $ and A 
respectively. For U < the stable modes of $ and A are at (ferromagnetism) and vr 
(superconductivity/CDW). Taking U > and expanding in q about the stable momenta 
gives the LG form: 






)A(g)-A(-g) 



(26) 



+ u^ Y. i^i • ^2) ($3 • '^-1-2-3) 

<li,<12,q3 

+ Mo E (Ai ■ A2) (A3 ■ A_i_2-3) 
91,92,93 

+ u, Y. (Ai ■ A2) ($3 ■ $-1-2-3) , 



91,92,93 

where $j = $(gi) (similarly for A) and: 
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Ud\{ix) 



3 9g^ 



U^ 



72 
216 

f/2 



E 



/'(e) 



^0 = 
1 

mo 
/(e) 



i + fx(o), 
f/5'x(o) 



+ 



3 dq^ 
1 ■ 
2? 



N{el 



Hn^Nie) 



E^iV(.) 
72 ^ e ^ ^ 



(27) 
(28) 

(29) 

(30) 

(31) 



where A^(e) is the non-interacting density of states. For U < the roles of $ and A are 
simply reversed. 

Equation E6| is the starting point for a renormalisation group (RG) analysis. We employ 



a one loop technique d la Wilson-Kogut [19,20|; i.e. successive elimination of 'fast-modes 



(large q) followed by rescaling of fields and momenta. This process renormalises the param- 
eters r^r^o and Mtt.cx- The analysis shows that the cross term, u^, is irrelevant, i.e. it flows 
to zero under the action of the RG. Therefore at the fixed point the spin and pseudospin 
are decoupled from each other. This implies that for U > the fixed point is the usual 
antiferromagnetic one, whilst it is S-CDW for U < 0. The sole effect of Ux is to renormalise 
Tc. Fig. |l] shows RG flows for different initial Ux but the same temperature, Tc{ux = 0). For 
Ux > the flow goes to the high temperature phase implying Tdux) < TciQi). Conversely 
for Ux < Q the ordered phase (AFM or S-CDW) is reached implying Tciux) > Tc{0). To 
conclude, Tc is depressed for positive Ux but increased for negative Ux- 

In summary, we have discovered a discrete group of symmetries and transformations of 
the Hubbard model on a bipartite lattice at half filling and enumerated its representations. 
Using these we have constructed a unique Hubbard Stratonovich decomposition of the Hub- 
bard interaction term, which treats spin and pseudospin on an equal footing. At a mean 
field level it is superior to other decompositions because it preserves all the transformation 
properties of the full theory. In addition, the unique decomposition we have derived has al- 
lowed us to construct a Landau- Ginzburg functional to examine the interplay of magnetism 
and superconductivity/CDW. For f/ > 0, the critical temperature for the onset of antifer- 
romagnetism, T^^^ , is found to be renormalised by ?7-paired fluctuations. Reciprocally, for 
U < 0, T^-*^^^ is renormalised by ferromagnetic fluctuations. Critical exponents are not 
affected. 
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FIGURES 
FIG. 1. RG flows for different initial values of Ux- For Ux = the full lines show flows starting 

infinitessinially above and below Tc- The broken lines are for the same temperature but non-zero 

initial Ux- 
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TABLES 



TABLE L Character Table for the group of C, f and Z. 
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